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ABSTRACT 


Six  different  methods  for  averaging  historical  salinity 
values  were  studied  for  their  applicability  in  the  determin¬ 
ation  of  depth  corrections  due  to  sound  speed  variations. 

The  averaging  techniques  range  in  complexity  from  the  use  of 
a  single  value  of  salinity  throughout  the  entire  water  column 
to  the  use  of  a  historical  salinity  average  profile  corrected 
for  the  surface  salinity.  Results  show  that  historically 
determined  salinity  can  be  used  in  the  computation  of  the 
depth  corrections  without  exceeding  the  accuracy  limits.  Two 
single  values  of  salinity  were  found  to  be  sufficient  to  cover 
the  West  Coast  of  the  United  States:  31  °/oo ,  applicable  north 
of  40°N,  and  33  %o ,  for  use  south  of  that  latitude. 
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I.  INTRODUCTION 


A.  THE  ECHOSOUNDER 

Despite  the  development  of  laser-  and  photobathymetry, 
the  echosounder,  or  fathometer,  is  still  the  most  commonly 
used  method  for  the  determination  of  depth  in  hydrography. 

The  transducer  emits  a  sound  pulse  which  is  reflected  by 
the  bottom  and  received  back  at  the  transducer.  The  time 
of  travel  of  this  sound  pulse  is  divided  by  two  and  this 
value  multiplied  by  the  assumed  speed  of  sound  in  sea- water,'*' 
thus  giving  the  depth  according  to  the  expression  z  =  v  t 
(distance  =  speed  x  time) .  This  transformation  is  made 
either  mechanically  or  electronically  within  the  fathometer 
itself,  and  the  result  displayed  is  the  nominal  or  fatho¬ 
meter  depth  beneath  the  transducer.  This  depth  is  not  equal 
to  the  true  depth,  since  the  assumed  sound  speed  generally 
does  not  equal  the  true  speed  of  sound  in  sea-water,  which 
varies  throughout  the  water  column.  One  example  of  an 
observed  sound  speed  profile  is  presented  in  Figure  1. 

For  determining  the  true  depth,  one  must  know  the  sound  speed 
profile  throughout  the  water  column  and  apply  it  according  to 
the  following  expression  [Greenberg  and  Sweers,  1972]  : 


Although  sound  speed  is  the  correct  term,  most  hydro- 
graphic  literature  uses  sound  velocity  instead. 
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In  this  equation,  v(zt)  represents  the  sound  speed  at  the  .  j 

depth  zt  where  the  signal  passes  at  the  time  t.  \ 

Some  echosounders  have  the  possibility  of  being  set  to  '  ] 

different  sound  speeds,  in  order  to  represent  as  closely  as 
possible  the  average  speed  of  sound  in  the  water  column.  In 

'  * 

most  cases,  however,  they  are  constructed  with  a  calibrated  f’ 

(■"•i 

speed  of  sound  assumed  in  the  device,  usually  1463.04  or  !- 

1500  m/s  (800  or  820  fathoms  per  second).  The  echosounders  j i 

used  by  the  National  Ocean  Survey  (NOS)  are  calibrated  for 

an  assumed  speed  of  sound  of  1463.04  m/s  which,  according  to  (• 

t 

Umbach  [1976]  is  a  value  reasonably  close  to  the  sound  speed 

in  most  waters  surveyed.  Corrections  are  then  applied  to  ; 

the  depths  obtained  by  the  fathometer  in  order  to  compensate  ,  1 

for  the  use  of  a  sound  speed  other  than  the  one  that  actually  i 

2 

affects  the  travel  of  the  sound  pulse. 

B.  ACCURACY  REQUIREMENTS  Fj 

The  International  Hydrographic  Bureau  (IHB)  established 
in  1968  the  following  requirements  for  accuracy  of  depth  !. 

i 

determination  [International  Hydrographic  Bureau,  1968]:  j 


Depth  Range 


Accuracy  Requirement 


0-11  fathoms  (0-  20  m)  1.0  ft 

11-55  fathoms  (20-100  m)  3.0  ft 

deeper  than  55  fathoms  1 %  of  the  depth 

The  above  standards  represent  the  maximum  allowable  error  in 
the  measurement  of  depth  due  to  all  sources,  including  the 
variations  of  sound  speed  within  the  water  column.  When  con¬ 
sidering  errors  due  to  sound  speed  alone,  no  errors  can  be 

equal  to  or  exceed  0.25%  of  the  depth,  which  means  that  the 

mean  sound  speed  must  be  known  to  within  ±4  m/s  [Umbach,  1976]. 
The  standards  mentioned  above  apply  only  to  conventional  echo¬ 
sounding  systems,  in  which  the  -sound  pulse  beam  is  vertical. 
Other  systems,  applying  oblique  acoustic  paths,  require  higher 
accuracy  standards. 

C.  SPEED  OF  SOUND  IN  THE  SEA 

Speed  of  sound  in  a  fluid  depends  on  the  density  (p)  ,  and 
compressibility  (k)  of  the  fluid  and  the  ratio  between  the 
specific  heats  of  the  fluid  at  constant  pressure  and  at  con¬ 
stant  volume  (y) ,  according  to  the  following  expression 
[Sverdrup,  Johnson  and  Fleming,  1942;  Bowditch,  1972]: 

v  = 

As  these  parameters  depend  on  temperature,  salinity  and 
pressure  (which  depends  on  depth) ,  it  is  more  convenient 


to  measure  these  properties  instead,  and  compute  the  sound 
speed  from  them. 


Different  expressions  have  been  developed  empirically 
by  Kuwahara  [1939],  Del  Grosso  [1950],  Wilson  [1960a,  b], 

Leroy  [1969]  and,  more  recently,  Clay  and  Medwin  [1977]. 
Mackenzie  [1960]  compared  several  of  these  formulas,  which 
showed  differences  between  them  of  less  than  3  m/s,  or  0.2% 
at  zero  depth  [Urick,  1975].  Speeds  used  for  depth  correc¬ 
tions  in  NOS  are  calculated  using  Wilson's  equation  [1960b], 
which  is  described  in  Appendix  A.  This  equation  is  the  most 
commonly  accepted  formula  due  to  its  large  range  of  applica¬ 
bility  and  its  relative  accuracy  (0.30  m/s  standard  deviation) 
[Wilson,  1960b,  1962]. 

D.  TRADITIONAL  METHODS  FOR  DETERMINING  DEPTH  CORRECTIONS 


I 


Several  different  procedures  can  be  used  to  determine 
the  depth  corrections,  or  correctors,  to  be  applied  to  echo- 
soundings.  These  corrections  can  either  be  found  directly 
or  by  knowing  the  sound  speed  profile.  The  sound  speed  can 
either  be  measured  directly  (direct  comparison  systems)  or 
computed  after  the  determination  of  some  characteristics  of 
the  water  column,  as  stated  above  (oceanographic  and  limno- 
logic  determination).  Some  of  the  different  procedures  are 
described  below: 


at  which  the  bar  is  lowered.  Because  these  are  comparative 
methods,  not  only  are  sound  speed  variations  included  in  the 
results  but  also  instrument  errors  [Umbach,  1976].  Both 
methods  are  limited  to  shallow  depths  and  good  weather  con¬ 
ditions  . 

a.  Leadline  Comparison  (Vertical  Cast) 

The  leadline  is  lowered  to  the  bottom  at  several 
depths  over  a  range  which  should  be  at  least  as  great  as  that 
likely  to  be  encountered  in  the  survey  area  [The  Hydrographer 
of  the  Navy,  1969].  Several  readings  are  taken  from  the 
fathometer  at  each  depth,  and  all  are  plotted  individually 
in  order  to  eliminate  random  errors.  In  areas  where  launch 
hydrography  joins  or  overlaps  ship  hydrography,  it  is  of 
particular  importance  to  use  this  method  on  both  ships  and 
launches  in  order  to  resolve  possible  discrepancies  in  the 
soundings  [Umbach,  1976].  This  method  is  limited  to  calm 
seas,  shallow  depths,  and  areas  of  hard  bottom  for  accurate 
readings . 

b.  Bar  Check 

The  depth  of  a  bar  indicated  by  the  echosounder 
is  compared  to  the  known  depth  at  which  the  bar  is  lowered 
below  the  transducer.  The  procedure  is  used  at  several  depths 
and  the  readings  are  taken  both  when  the  bar  is  lowered  and 
when  it  is  raised.  Since  the  bar  is  handled  over  the  side, 
this  method  is  usually  limited  to  small  vessels  and  shallow 
depths,  although  it  is  sometimes  performed  mechanically  on 
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safer: 


large  vessels.  The  depth  limitation  still  holds,  however, 
due  to  difficulties  in  maintaining  vertical  lowering  lines 
and  in  keeping  the  bar  within  the  beam  of  the  transducer. 
This  procedure  should  only  be  used  in  calm  seas  with  light 
winds  and  no  strong  current. 

2 .  Direct  Measurement  of  Sound  Speed 

Sound  speed  can  be  measured  directly  by  lowering  a 
sound  velocimeter  to  various  depths  and  recording  the  speeds 
indicated.  Although  accuracies  obtained  by  this  method 
can  be  of  the  order  of  a  few  tenths  of  a  meter  per  second 
[Ulonska,  1972],  sound  velocimeters  are  not  generally  used 
because  they  are  very  expensive  instruments  [Calder,  1975]. 

3.  Measurements  of  Salinity,  Temperature  and  Depth 


These  parameters  can  be  determined  by  direct  meas¬ 
urement  either  with  Nansen  casts  (temperature  and  salinity 
at  specific  depths)  or  with  continuous  sampling  instruments 
which  measure  salinity  (or  conductivity) ,  temperature  and 
depth  (STD  or  CTD) .  The  latter  instruments  can  also  be  of 
an  expendable  type.  Sound  speed  at  various  depths  is  then 
computed  using  an  equation  relating  salinity,  temperature 
and  depth  (pressure)  to  sound  speed,  such  as  Wilson's  equa¬ 
tion  [1960b],  or  they  can  be  extracted  from  tables,  wherein 
sound  speed  is  tabulated  for  different  values  of  temperature, 
pressure  and  salinity  [Heck  and  Service,  1924;  Matthews, 

1927,  1939;  Kuwahara,  1939]. 
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4.  Historical  Tables  or  Atlases 


Sound  speed  or  salinity  and  temperature  data  have 
been  tabulated  for  various  areas  and  seasons.  Matthews' 
tables  [1939]  provide  corrections  to  be  applied  to  the  depths 
obtained  with  fathometers  calibrated  to  a  sound  speed  of  1463 
or  1500  m/s.  These  depth  corrections  are  given  as  a  function 
of  nominal  depth  and  are  tabulated  for  52  different  oceanic 
areas  identified  by  Matthews  [Greenberg  and  Sweers,  1972]. 

This  method  is  not  used  by  NOS  but  it  is  used  in  some  coun¬ 
tries  for  depths  exceeding  200  m.  For  example,  this  is  the 
only  method  which  Australia  uses  for  determining  corrections 
in  water  depths  of  more  than  200  m  [Compton,  1980]. 

E.  THE  LAYER  METHOD 

1 .  Options 

In  addition  to  the  above  techniques  to  establish  the 
depth  correction,  yet  another  more  sophisticated  method  is  in 
widespread  use--the  "summation  of  layers"  method  [Umbach, 

1976],  This  method  is  applicable  when  the  sound  speed  is 
determined  either  by  direct  measurement  or  by  computation  from 
the  depth,  temperature  and  salinity  parameters  of  the  water 
column.  This  method  can  be  applied  in  two  ways,  which  are 
described  in  the  "Hydroplot/Hydrolog  Systems  Manual"^  [Wallace, 
1971]  : 

^  The  Hydroplot/Hydrolog  system  is  an  automated  data  aqui- 
sition  and  processing  system  in  use  by  NOS,  respectively  with 
and  without  a  plotting  capability.  A  description  of  this 
system  is  given  in  the  Hydrographic  Manual,  Fourth  Ed.  [Umbach, 
1976]. 


i- 


18 


a.  Curve- fit  Method 

The  water  column  is  partitioned  into  layers  of 
varying  thickness  which  are  predetermined  by  the  operator. 
Layer  mid-depths  are  then  computed  and  the  sound  speed  at 
those  depths  interpolated  either  graphically  or  mathemati¬ 
cally  from  a  curve  passed  through  the  points  at  which  the 
sound  speed  has  been  computed. 

b.  No-curve-fit  Method 

The  water  column  is  partitioned  as  above  but  the 
depths  at  which  the  observations  have  been  made  are  consid¬ 
ered  as  the  mid-depths,  the  boundaries  between  layers  being 
defined  as  half-way  between  two  successive  observations  of 
the  cast.  One  example  illustrating  the  no-curve-fit  method 
is  shown  in  Figure  2. 

2 .  Computation  of  the  Depth  Corrections 

a.  Layer  Sound  Speed 

For  each  layer  a  sound  speed  is  determined  for 
the  corresponding  mid-depth  by  any  of  the  methods  described 
above.  This  value  of  sound  speed  is  assumed  to  apply  through¬ 
out  the  entire  layer. 

b.  Layer  Factor 

A  factor  is  computed  for  each  layer,  based  on  the 
sound  speed  for  which  the  instrument  is  calibrated,  using  the 
following  expression  [Umbach,  1976]: 


f actor ^ 
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where  is  the  actual  sound  speed  for  the  i-th  layer  which 
has  been  determined  for  the  corresponding  mid-depth  (see 
preceding  paragraph) ,  and  C  is  the  assumed  sound  speed  of 
the  instrument. 

c.  Layer  Correction 

The  above  factor  is  multiplied  by  the  layer 
thickness,  Az^  (the  difference  between  the  depth  at  the 
bottom  and  at  the  top  of  the  layer)  in  order  to  obtain  the 
correction  to  the  depth  (LC^)  due  to  the  variation  in  sound 
speed  of  that  layer: 

LC ■  =  factor,  x  Az. 

1  11 

The  water  column  above  the  transducer  is  eliminated  when 
determining  the  corrections  for  the  upper  layer. 

d.  Bottom-of- layer  Correction 

The  layer  corrections  are  added  from  the  surface 
downward  in  order  to  obtain  the  total  correction  at  the  bottom 
of  each  layer: 


i 


where  is  the  correction  at  the  bottom  of  the  i-th  layer, 

e.  Fathometer  Depth  Correction  Table 

The  above  corrections  correspond  to  the  true 
depth  of  the  bottom  of  each  layer.  To  obtain  the  fathometer 
depth  to  which  these  corrections  will  apply,  one  must  sub- 
the  correction  (C^)  from  the  true  depth  of  the  bottom 


m. : 
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tract 


of  the  corresponding  layer  (z^  : 


where  Cp^  is  the  fathometer  depth  corresponding  to  the  i-th 
1 

] ayer . 

f.  Correction  Curve 

The  values  of  the  corrections  are  plotted  against 
the  fathometer  depth.  A  curve  is  drawn  to  fit  these  depth- 
correction  pairs.  From  this  curve  the  depth  correction  appli¬ 
cable  to  any  reading  of  the  fathometer  (within  the  range  for 
which  the  curve  has  been  computed)  can  be  obtained.  An  example 
of  such  a  curve,  for  which  a  transducer  at  lero  depth  has  been 
considered,  is  shown  in  Figure  5. 

F.  VARIATION  OF  THE  PARAMETERS 

The  variation  of  sound  speed  with  depth,  temperature 
and  salinity  depends  on  the  initial  conditions  of  these  par¬ 
ameters.  For  a  temperature  of  14°C,  salinity  of  34  °/oo  and  at 
zero  depth  (sound  speed  under  these  conditions  is  1502.94  m/s) 
an  increase  of  one  degree  in  temperature  will  increase  sound 
speed  by  about  4.5  m/s;  an  increase  of  1  %0  in  salinity  will 
increase  the  sound  speed  by  about  1.3  m/s;  an  increase  in 
depth  of  100  m  (about  10  atmospheres  in  pressure)  will  in¬ 
crease  sound  speed  by  about  1.7  m/s  [Urick,  1975].  Latitude 
also  has  some  influence,  since  pressure  increases  from  the 
equator  to  the  poles  because  of  the  increase  of  the  gravita¬ 
tional  attraction.  This  influence,  however,  is  very  small 
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and  can  be  ignored.  An  example  of  the  variation  of  sound 
speed  with  temperature,  salinity  and  depth  is  shown  in 
Figure  1. 

Temperature  is  the  parameter  that  has  the  largest  influ¬ 
ence  on  sound  speed.  Also,  it  shows  more  variability  with 
time.  (Jmbach  [1976]  states  that,  for  the  sound  speed  to  be 
accurate  to  within  ±4  m/s,  temperature  and  salinity  must  be 
determined  with  an  accuracy  of  +1°C  and  1  °/oo  ,  respectively. 
However,  recent  studies  have  been  performed  by  the  Testing 
Division,  Office  of  Marine  Technology  (OMT),  NOAA,  regarding 
the  influence  of  the  different  parameters  on  sound  speed. 

These  studies  show  that  if  temperatures  can  be  determined 
with  an  accuracy  of  ±.1°C,  the  salinity  accuracy  need  only 
be  within  ±.  3  °/oo  to  still  meet  the  requisite  accuracy  standards 
[Bivins,  1976]. 

G.  USE  OF  HISTORICAL  INFORMATION 

Historical  information  has  been  used  for  some  time  in  the 
computation  of  depth  corrections  for  echosoundings ,  but  usually 
only  for  depths  in  excess  of  200  m  [Sherwood,  1974;  Urick, 
1975].  Yeager  [1979]  studied  whether  depth  corrections  of 
sufficient  accuracy  could  be  determined  from  historical  salin¬ 
ity  and  temperature  only.  He  chose  an  area  which  is  repre¬ 
sentative  of  the  near  shore  region  along  the  East  Coast  of 
the  United  States.  The  historical  sound  speed  data  there 
indicated  that  the  temporal  and  spatial  variability  of  tem¬ 
perature  exceeded  the  acceptable  limits  for  sounding 
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corrections,  thus  precluding  the  use  of  historical  data  alone 
for  correction  determination.  In  his  study  he  also  concluded 
that,  in  that  area  (a  little  more  than  one  degree  square  with 
depths  less  than  200  m)  salinity  values  were  sufficiently 
stable  to  allow  acceptable  depth  corrections  to  be  computed 
from  historically  derived  salinity  values  combined  with  in 
situ  temperature  measurements. 

In  Australia  data  have  been  collected  over  the  contin¬ 
ental  shelf  since  1965  which  show  salinity  to  vary  generally 
between  34.5  and  36  %o ,  changing  only  slightly  with  depth 
down  to  500  m  [Calder,  1975J.  At  present,  the  Australian 
Hydrographic  Service  exclusively  uses  a  value  of  35  °/oo  for 
the  calculation  of  sound  speed  in  depths  less  than  200  m 
[Compton,  1980] . 

H.  OBJECTIVES 

The  purpose  of  this  thesis  is  to  analyze  the  possibility 
of  using  historical  salinity  data  averages  for  a  large  area 
with  in  situ  temperature  measurements  to  compute  depth  cor¬ 
rections  and  still  remain  within  the  required  accuracy  limits. 
The  historical  salinity  data  were  treated  using  several  dif¬ 
ferent  techniques  which  increased  in  accuracy  as  more  complex 
averaging  schemes  were  used. 

Temperature-salinity  observations  generally  have  to  be 
made  at  least  once  a  month  and  be  of  sufficient  number  to 
cover  the  entire  area  sounded  [Umbach,  1976].  The  use  of 
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historical  salinity  data  and  XBT  (expendable  bathythermo¬ 


graph)  temperature  information  could  save  a  significant 
amount  of  time  and  money  if  results  remain  within  the 
accuracy  standards.  Also,  since  little  time  is  needed  for 
each  XBT  observation,  more  frequent  temperature  measurements 
could  be  made  (both  in  time  and  space),  thus  increasing  the 
accuracy  of  the  corrections  used  throughout  the  survey. 


II.  PROCEDURE  AND  METHODOLOGY 


Some  of  the  terms  to  be  used  in  the  following  chapters 
are  defined  below: 

Historical  data:  data  resultant  from  some  averaging 
scheme  applied  to  previous  observations  in  a  large 
area . 

Inferred  data:  data  obtained  from  both  historical  and 
observed  salinities,  used  to  represent  the  data  for  a 
large  region.  It  has  not  been  measured,  and  consists 
of  a  single  value. 

Station:  group  of  measurements  taken  at  a  particular 

time  in  a  particular  location.  Each  station  consists 
of  several  observations  taken  over  a  given  depth 
interval . 

Observation:  measurement  taken  at  a  particular  depth 

in  one  station. 

True  correction:  depth  correction  resultant  from  the 
application  of  the  observed  data. 

Estimated  correction:  depth  correction  for  which  the 
computation  of  some  historical  data  were  used. 

Error:  difference  between  the  estimated  and  true 

corrections . 

Percentage  of  error:  division  of  the  error  by  the 
depth  over  which  the  correction  applies. 


A.  STUDY  AREA 

The  area  selected  for  this  study  was  the  continental  shelf 
region  of  the  West  Coast  of  the  United  States,  extending  off¬ 
shore  generally  to  depths  of  about  1500  m.  The  study  area 
can  be  seen  in  Figure  4.  This  area  was  chosen  on  the  basis 
of  several  factors: 
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Figure  4.  Study  area,  showing  the  inshore  limits  of  the 
Alaskan  and  Californian  water  masses. 
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1.  It  represents  a  fairly  large  area  wherein  a  consider¬ 
able  amount  of  data  were  available.  For  the  purpose  of  this 


study,  a  large  area  was  desired  to  generalize  the  applica¬ 
bility  of  the  results. 

2.  Most  hydrographic  surveys  are  conducted  over  the 
continental  shelf  region. 

3.  The  area  was  well  suited  for  acquisition  of  histor¬ 
ical  data,  as  will  be  seen  later  in  this  chapter. 

4.  The  extension  to  1500  m  depths  guaranteed  that  all 
shelf  features,  including  submarine  canyons,  would  be  con¬ 
sidered  . 

The  study  area  was  subdivided  into  three  sub-areas: 
Northern,  Southern  and  Strait  of  Juan  De  Fuca  and  Puget  Sound. 
This  latter  area  was  considered  separately  because  of  its 
particular  salinity  characteristics.  The  North-South  division 
of  the  Pacific  Ocean  waters  along  latitude  40°N  corresponds 
to  the  separation  of  the  two  different  water  masses  consid¬ 
ered  in  the  historical  salinity  file. 

B.  DATA  SOURCES  AND  TREATMENT  OF  DATA 

1 .  Historical  Data 

Historical  data  for  this  study  were  obtained  through 
the  PROFIT  subroutine  of  the  Integrated  Command  ASW  Predic¬ 
tion  System  (ICAPS).  The  ICAPS  water  mass  history  file  was 
developed  by  the  U.S.  Naval  Oceanographic  Office  (NAVOCEANO) 
and  contains  historical  information  on  the  water  masses  of 
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the  Northern  Hemisphere  and  Indian  Ocean  [Department  of  the 
Navy,  1978].  Each  ocean  is  divided  into  several  areas  and 
subareas.  In  each  of  these  subareas  up  to  five  different 
water  masses  may  be  included.  Figure  5  shows  the  North 
Pacific  Ocean,  Pacific  Area  A  divisions  and  the  water  masses 
present  in  Area  A.  The  water  mass  files  are  also  grouped 
seasonally:  winter  (January  through  March) ,  spring  (April 

through  June),  summer  (July  through  September),  fall  (October 
through  December) . 

Within  the  study  area  the  ICAPS  water  mass  file  con¬ 
sists  of  only  two  water  masses:  Alaskan  water  (north  of  40°N 
and  Californian  water  (south  of  40°N) .  These  are  shown  in 
Figure  5(b).  The  inshore  limits  of  these  water  masses  are 
shown  in  Figure  4.  The  Alaskan  water  mass  was  used  for  the 
near  shore  area  off  Washington,  although  this  area  is  not 
included  in  either  of  the  two  water  masses. 

For  ICAPS  data  retrieval  and  processing,  several  dif¬ 
ferent  program  modules  are  available  [Chace,  1979].  In  this 
study  only  the  PROFIL  subroutine  was  used.  One  of  the  out¬ 
puts  of  PROFIL  is  historical  salinity.  It  was  used  to  obtain 
information  for  each  area  during  all  four  seasons.  The 
salinity  profiles  from  these  eight  runs  were  punched  onto 
cards  and  are  shown  in  Table  I. 

2 .  Observed  Data 

The  main  source  of  observed  data  was  the  Fleet  Numer¬ 
ical  Oceanography  Center  (FNOC)  file  entitled  MOODS  (Master 
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Figure  5(a) 


Figure  5(b) 
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Figure  5(a)  North  Pacific  Ocean  locator  chart,  (b)  Pacific 
area  A,  and  (c)  Pacific  area  A  water  masses,  as  defined 
in  the  ICAPS  history  file. 
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HISTORICAL  SALINITY  FOR  THE  WEST  COAST 
OF  THE  UNITED  STATES,  FROM  ICAPS 
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Oceanographic  Observation  Data  Set).  Some  additional  data 
were  obtained  from  the  Pacific  Marine  Center  (PMC),  NOAA, 
Seattle,  Washington.  The  MOODS  file  contains  data  from  sea 
water  observations  all  over  the  world.  Different  types  of 
observations  are  included:  STD,  CTD,  Nansen  cast,  BT  or  XBT, 
and  others.  Each  station  contains  the  geographical  location, 
date,  time  and  a  code  for  the  ship/agency  that  made  the  ob¬ 
servation.  Temperature,  salinity,  sound  speed,  oxygen  con¬ 
centration  or  other  properties  are  also  provided  depending 
on  the  type  of  observation.  The  depth  interval  between 
observations  for  different  stations  is  not  standardized. 

The  time  interval  chosen  for  the  study  was  the  ten 
year  period, from  1964  to  1973.  A  more  recent  period  was  not 
chosen  because  a  very  large  percentage  of  the  recent  data  is 
from  XBT  observations,  which  does  not  provide  salinity  infor¬ 
mation. 

A  direct  extract  from  the  MOODS  file  was  not  readily 
usable  on  the  Naval  Postgraduate  School  computer  since  the 
MOODS  data  are  stored  with  a  variable  block  size.  A  nine 
track  tape,  1600  bits  per  inch,  constant  block  size,  was  then 
prepared  at  FNOC  which  contained  temperature  and  salinity  for 
depths  standardized  to  0,  30,  60,  90,  120,  150,  175,  200,  250 
and  300  m.  Values  in  excess  of  300  m  were  not  given,  even  if 
they  had  existed  in  the  original  data.  In  the  meantime,  some 
data  were  punched  onto  cards,  extracted  from  a  printout  ob¬ 
tained  from  the  original  data.  Since  this  data  set  had  a 


better  profile  definition,  especially  useful  for  the  inter¬ 
mediate  values  between  the  surface  and  50  m,  it  was  decided 
to  use  this  punched  data  in  substitution  for  the  correspond¬ 
ing  stations  on  the  tape. 

A  total  of  3459  stations  were  analyzed,  with  the 
distribution  shown  in  Table  II.  For  each  of  these  stations 
a  depth  correction  was  calculated  based  on  the  procedures 
outlined  in  the  previous  chapter.  Such  corrections  are  con¬ 
sidered  true  depth  corrections  and  are  the  values  against 
which  the  estimated  depth  corrections  (based  on  some  mean 
historical  salinity  value)  must  be  compared. 

TABLE  II 

DISTRIBUTION  OF  THE  STATIONS  BY  AREA  AND  SEASON 


Winter 

Spring 

Summer 

Fall 

Total 

North 

248 

679 

643 

195 

1765 

South 

450 

350 

229 

195 

1222 

Puget  Sound/ 
Strait  of 

127 

149 

100 

96 

472 

Juan  De  Fuca  - 

Total  3459 


Appendix  B  shows  all  the  station  locations  for  each 
area  and  season.  Multiple  occupations  of  a  station  may  have 
been  made  which  are  not  indicated  by  these  figures. 
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C.  DATA  ANALYSIS 


1 .  Salinity 

From  the  ICAPS  data  file,  mean  salinity  values  were 
computed  from  the  surface  to  a  given  depth,  for  each  area 
and  season.  This  mean  was  determined  not  only  for  the  depth 
of  each  given  value  of  historical  salinity  (Table  I)  but  also 
for  the  mid-depths  between  two  successive  values.  Table  III 
contains  the  values  obtained.  These  data  were  then  used  to 
establish  a  variety  of  historical  salinity  averages  to  be 
described  below. 

Estimated  depth  corrections  were  then  determined  for 
each  station  using  the  observed  temperature  and  an  historical 
salinity  value  derived  from  one  of  four  different  methods: 

a.  Gross  Historical  Salinity  Average  (Method  1) 

This  method  uses  a  single  value  of  salinity 

throughout  the  entire  water  column.  This  salinity  value 
was  determined  by  averaging  the  historical  salinity  data 
from  the  surface  to  200,  1000  or  3000  m.  The  choice  of  the 
depth  depended  on  whether  the  deepest  depth  for  the  station 
was  less  than  200,  between  200  m  and  3000  m,  or  greater  than 
3000  m. 

b.  Fine  Historical  Salinity  Average  (Method  2) 

This  method  was  similar  to  the  previous  one 

except  that,  instead  of  three  possible  values  for  salinity, 
ten  possible  values  were  considered:  averages  from  the  sur¬ 
face  to  50,  100,  150,  200,  300,  500,  750,  1000,  1500  and 
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TABLE  III 


HISTORICAL  SALINITY  AVERAGES  FOR  THE 
WEST  COAST  OF  THE  UNITED  STATES 
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3000  m  depth.  As  before,  the  depth  greater  than  or  equal  to 
the  deepest  depth  observed  was  chosen.  For  example,  if  the 
depth  were  250  m,  the  salinity  average  would  be  the  one  cor¬ 
responding  to  the  300  m  water  column. 

c.  Historical  Salinity  Curve  (Method  3) 

The  value  from  the  historical  salinity  curve, 
which  corresponded  to  a  given  depth,  was  obtained  by  interpo¬ 
lation  from  the  existing  values  extracted  from  the  ICAPS  file. 

d.  Adjusted  Historical  Salinity  Curve  (Method  4) 

Since  the  salinity  value  observed  at  the  sea  sur¬ 
face  can  sometimes  deviate  considerably  from  the  historical 
surface  value,  an  alternative  to  the  historical  salinity  curve 
was  considered.  A  linear  adjustment  was  made  to  this  curve 

so  that  the  surface  value  would  equal  the  observed  salinity 
value,  and  the  value  at  30  m  would  equal  the  historical  data 
(ICAPS)  at  that  depth.  Below  30  m  the  historical  curve  was 
used. 

Salinity  values  illustrating  these  first  four  methods 
as  applied  to  a  particular  station  are  shown  in  Figure  6.  In 
light  of  the  results  obtained  using  these  four  methods,  two 
new  approaches  were  taken: 

e.  Inferred  Value  (Method  5) 

A  single  value  of  salinity  was  used  for  each  area, 
independent  of  season  and  depth.  This  value,  which  was  inferred 
from  both  the  historical  data  and  the  actual  observations,  was 
judged  to  represent  average  conditions  fpr  the  area. 
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f.  Inferred  Value  With  Observed  Surface  Salinity 
(Method  6) 

This  method  is  similar  to  the  previous  method, 
except  that  the  salinity  which  is  applied  to  the  first  layer 
is  the  observed  surface  salinity.  This  method  is  also  simi¬ 
lar  to  method  4,  except  it  only  uses  one  value  of  salinity 
throughout  the  column  for  depths  other  than  the  surface. 

2 .  Corrections 

Each  of  the  six  methods  was  applied  as  in  a  normal 
survey.  Estimated  depth  corrections  were  computed  for  each 
layer  using  the  no-curve-fit  method,  i.e.,  the  depth  of  the 
observed  temperature  was  used  as  the  mid-depth  for  the  layer. 
The  procedures  indicated  in  section  I.E.  were  followed  using 
steps  a  to  d; 

(aj  Estimated  sound  speed  for  each  layer. 

(b)  Computation  of  the  corresponding  factor. 

(c)  Layer  correction. 

(d)  Bottom  of  layer  correction 

The  correction  curves  for  the  observed  temperature  and  salin¬ 
ity  values  and  those  based  upon  each  of  the  first  four  methods 
are  shown  in  Figure  7.  The  estimated  correction  at  the  bottom 
of  each  layer  was  compared  to  the  "true"  correction  wnich  was 
obtained  using  the  observed  salinity. 

Error.  =  estimated  correction.  -  true  correction. 

1  i  l 

The  differences  were  divided  by  the  depth,  in  order  to  obtain 
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Correction  curves  for  the  first  four  methods. 


the  percentage  of  error  for  that  particular  depth: 

Error.. 

Percentage  of  error,  =  - - 

1  zi 

where  is  the  depth  of  the  bottom  of  the  layer.  If  this 
value  were  greater  than  or  equal  to  0.251,  the  method  was 
considered  unsuitable  at  that  depth,  since  the  accuracy  re¬ 
quirements  could  not  be  met.  The  values  compareu  were  the 
corrections  at  the  true  depth  at  the  bottom  of  each  layer. 

In  reality,  the  values  at  the  fathometer  depth  (true  depth 
minus  correction)  should  have  been  compared.  This  was  not 
done  since  the  fathometer  depths  would  probably  be  different 
for  each  error.  However,  the  resultant  error  from  this  approx¬ 
imation  is  insignificant- -of  the  order  of  a  few  millimeters 
(see  Figure  8) . 

All  the  computations  for  the  first  four  methods  were 
executed  sequentially  in  the  same  program,  shown  in  Appendix 
C.  This  program  was  also  used  for  the  fifth  and  sixth  methods. 
An  output  of  this  program  for  a  particular  station  is  shown  in 
Figure  9.  The  meaning  of  the  different  headings  is  as  follows: 
SHI:  gross  historical  salinity  average  (method  1) 

SH2:  fine  historical  salinity  average  (method  2) 

Z  :  depth  of  the  observation 

ID  :  depth  at  the  bottom  of  each  layer 

T  :  observed  temperature 

S  :  observed  salinity 
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figure  8(a)  Figure  8(b)  Figure  8(c) 

Figure  8.  Depth  corrections  as  a  function  of  "true  depth"  (a)  and 
"true  depth  minus  correction"  (b)  for  true  salinity  and  method 
Figure  8(c)  illustrates  the  difference  between  the  two  methods. 
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Figure  9.  Program  output  for  a  particular  station. 


HS  :  historical  salinity  value  (method  3) 

HS' :  adjusted  historical  salinity  curve  (method  4) 

V  :  sound  speed  computed  from  Z,  T  and  S 

Vl-4:  sound  speed  computed  using  methods  1-4 

COR  :  correction  applied  at  depth  ID  computed 
from  V 

C0R1-4:  correction  obtained  using  methods  1-4 

DCOR1-4:  error  of  methods  1-4  (difference  between 

C0R1-4  and  COR) 

PCI-4:  percentage  of  error  of  methods  1-4 

(DC0R1-4  divided  by  ID) 


i 
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III.  DISCUSSION  OF  THE  RESULTS 
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Based  on  the  foregoing  analysis,  four  different  possi¬ 
bilities  have  been  considered  for  each  of  the  six  salinity 
averaging  methods  tested: 

(1)  No  error  exceeds  the  limits  at  any  given  depth. 

(2)  The  error  limit  is  exceeded  only  in  the  surface 
layer . 

(3)  The  error  limit  is  exceeded  between  the  surface 
layer  and  the  layer  which  contains  the  30  m  depth. 

(4)  Excessive  errors  have  been  found  for  layers  deeper 
than  30  m. 

A.  USE  OF  HISTORICAL  SALINITY  DATA  DERIVED  FROM  ICAPS 

(FIRST  FOUR  METHODS) 

The  results  obtained  from  using  historical  salinity  data 
derived  from  the  ICAPS  water  mass  file  are  shown  in  Table  IV. 
A  discussion  of  these  results,  categorized  by  area,  follows. 

1 .  Northern  Area 

The  large  variation  of  surface  salinity  in  this 
region  caused  the  depth  correction  for  the  first  layer  to 
exceed  the  acceptable  limits  for  the  first  three  salinity 
averaging  methods  considered.  Excessive  errors  were  found 
for  31.6%,  28.0%  and  24.4%  of  the  stations  examined,  respec¬ 
tively.  A  slight  improvement  is  noted  between  methods  1  and 
3,  indicating  that  some  benefit  is  derived  from  the  use  of 
an  historical  mean  salinity  profile.  However,  when  the 
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TABLE  IV 


ERRORS  RESULTANT  FROM  THE  USE  OF  HISTORICAL  SALINITY 
DERIVED  FROM  ICAPS  (METHODS  1  TO  4). 

Results  are  shown  for  the  surface  layer  and  for  water 
depths  shallower  or  deeper  than  30  m.  Numbers  indicate 
the  number  of  stations  which  exceeded  the  accuracy 
standards.  Percentages  are  indicated  in  parenthesis. 
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historical  salinity  profile  was  adjusted  to  match  the  surface 
value  (method  4) ,  there  were  no  stations  which  exceeded  the 
acceptable  limits. 

The  depth  error  of  the  first  layer  was  often  so  large 
that  sometimes  the  offset  created  induced  an  unacceptably 
large  error  in  the  second  layer.  This  occurred  for  14.9%, 

11.6%  and  10.2%  of  the  stations  examined,  respectively.  In 
very  few  cases  (2.7%,  1.9%  and  1.5%)  the  excessive  difference 
in  surface  salinities  created  an  excessive  error  which  extended 
to  the  third  layer. 

There  are  two  main  reasons  why  this  error  is  so  large: 
1)  the  observed  salinity  extends  to  very  low  values,  and  2) 
most  of  the  data  only  are  defined  at  the  surface  and  50  m 
isopleths  which  creates  a  very  wide  first  layer  (15  m)  over 
which  that  value  is  applied.  The  errors  which  occur  due  to 
low  surface  salinities  generally  do  not  cause  the  accuracy 
limits  to  be  exceeded  except  close  to  the  surface,  for  depths 
rarely  exceeding  30  m.  The  salinity  may  be  very  low  at  the 
surface  but  it  increases  quite  rapidly  to  nominal  oceanic 
values  in  a  few  meters. 

The  error  due  to  low  surface  salinities  is  of  par¬ 
ticular  importance  during  the  spring  season  (April-June)  due 
to  the  large  outflow  of  the  rivers  of  Northern  California  and 
Oregon,  especially  the  Columbia  River.  It  is  also  noticeable 
in  the  winter  season,  but  is  considerably  less  important  in 
summer.  An  example  where  this  situation  frequently  occurred 


was  in  the  area  of  influence  of  the  Columbia  River,  where 
on  one  occasion  a  salinity  as  low  as  8.02  °/6o  was  observed 
during  the  spring  season. 

2 .  Southern  Area 

Considering  the  least  sophisticated  of  the  salinity 
averaging  schemes  (method  1),  only  five  stations  out  of  1222 
(0.4%)  showed  any  kind  of  deviation  which  exceeded  the  accept¬ 
able  depth  limits.  Of  these  five,  only  one  station,  located 
adjacent  to  the  Golden  Gate  Bridge  in  San  Francisco,  exceeded 
those  limits  in  a  layer  other  than  the  surface.  For  this 
station,  unacceptable  depth  errors  were  found  at  depths 
between  70  and  80  m  due  to  an  inversion  in  the  salinity  pro¬ 
file.  Of  the  other  four  station',  only  one  exceeded  the 
error  limit  when  the  second  method  was  used.  All  these  sta¬ 
tions  were  located  in  water  more  than  300  m  deep,  except  one 
located  about  8  miles  south  of  Point  Reyes,  about  half-way 
to  the  Farallon  Islands.  None  of  these  observations  exceeded 
the  limits  for  the  third  and  fourth  methods.  The  general 
pattern  resulted  in  a  slight  improvement  in  accuracy  from  the 
first  to  the  third  methods,  and  a  significant  improvement 
when  the  fourth  method  was  applied.  However,  there  were  cases 
where,  because  of  the  large  range  in  salinity,  method  1  was 
more  accurate  than  method  3.  The  pattern  of  the  accuracy  of 
the  four  methods  was  similar  to  that  experienced  in  the 
Northern  area. 
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3 .  Strait  of  Juan  De  Fuca  and  Puget  Sound 

Although  method  4  demonstrated  a  significant  increase 
in  accuracy,  none  of  the  four  methods  which  used  historical 
salinity  derived  from  the  ICAPS  file  was  deemed  accurate 
enough.  This  file  does  not  define  any  water  mass  applicable 
in  the  Strait  of  Juan  De  Fuca  and  Puget  Sound,  and  the  extrap¬ 
olation  of  the  ICAPS  data  to  fit  this  region  proved  to  be 
insufficient  in  accuracy. 

B.  USE  OF  INFERRED  SALINITY  VALUE 

(FIFTH  AND  SIXTH  METHODS) 

Salinity  values  of  31,  33  and  29  °/oo  were  selected  as 
representative  of  the  Northern,  Southern  and  Strait  of  Juan 
De  Fuca/Puget  Sound  areas,  respectively.  These  values  were 
applied  both  throughout  the  entire  water  column  (method  5) 
or  else  were  replaced  by  the  observed  surface  salinity  value 
for  the  first  layer  only  (method  6).  The  number  of  stations 
that  exceeded  the  acceptable  limits  for  each  of  these  methods 
is  shown  in  Tables  V  and  VI.  A  further  analysis  of  these 
results  follows. 

1 .  Northern  Area 

Of  the  176S  stations  examined,  323  stations  (18.3a) 
had  a  surface  layer  error  that  exceeded  the  acceptable  limits 
when  a  value  of  31  °/oo  was  used  throughout  the  entire  water 
column  (method  5).  Only  117  (6.61)  exceeded  the  limits  in 
deeper  layers,  and  only  one  in  a  layer  deeper  than  30  m. 


TABLE  V 


ERRORS  RESULTANT  FROM  THE  USE  OF  A  SINGLE  VALUE  OF  SALINITY 
THROUGHOUT  THE  ENTIRE  WATER  COLUMN  l METHOD  5)* 


*Table  descriptions  as  in  Table  IV. 


TABLE  VI 

ERRORS  RESULTANT  FROM  THE  USE  OF  A  SINGLE  VALUE  OF  SAL  IN  Hi 
THROUGHOUT  THE  ENTIRE  WATER  COLUMN  (METHOD  o)* 
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*Table  descriptions  as  in  Table  IV 


SO 


When  the  observed  surface  saiinitv  value  was  used  in  the 


first  laver,  no  stations  exhibited  an  unacceptable  error. 

2 .  Southern  Area 

No  excessive  errors  were  found  when  either  method  5 
or  o  was  applied,  for  which  a  constant  value  of  3  3  °/co  was 
assumed  for  this  area. 

3 .  Strait  of  Juan  Pe  Fuca  and  Puget  Sound 

A  salinity  value  of  19  Voo  was  assumed  for  this  area. 
Of  the  4 " 1  stations  examined,  113  stations  123.91)  exceeded 
the  limits  when  this  value  was  used  by  itself,  but  only  23 
.4.91)  when  it  was  substituted  in  the  first  layer  by  the 
observed  > u r  f a c e  s a  1  i n  1 1 v  value. 

It  ha.'  been  seen  that  for  the  three  areas,  and  for  any 
of  the  methods  considered,  most  of  the  errors  exceeding  the 
accuracy  Limits  occurred  in  t ho  sort  ice  laver  on  1  v.  This  does 
not  affect  the  accuracy  standards  when  the  depth  is  deeper 
than  about  30  m,  which  represents  most  of  the  area  over  which 
hydrographic  surveys  are  conducted. 


IV.  CONCLUSIONS 


The  results  shown  in  the  previous  chapter  indicate  that 
historical  data  can  indeed  be  used  in  the  computation  of 
depth  corrections  due  to  variations  in  sound  speed  for  the 
West  Coast  of  the  United  States.  Data  derived  from  the  ICAPS 
salinity  file  are  completely  acceptable  only  in  the  Southern 
area,  although  good  approximations  can  be  made  in  the  Northern 
area.  The  situation  in  the  Northern  area,  however,  can  be  an 
artifact  created  by  the  observed  data  format.  If  the  surface 
value  as  observed  in  the  sea  is  applied  to  adjust  the  histor¬ 
ical  salinity  curve  from  ICAPS,  then  the  results  are  completely 
correct  for  all  depths  and  seasons,  both  in  the  Northern  and 
Southern  areas.  This  by  itself  would  indicate  sufficient 
confidence  to  apply  this  method  to  actual  hydrographic  sur¬ 
veys  in  the  areas  considered.  It  would  only  require  the  use 
of  a  table,  such  as  Table  I,  and  the  measurement  of  the  sur¬ 
face  salinity,  the  temperature  being  measured  with  an  XBT. 

The  advantages  of  this  method  over  what  is  presently  being 
done  are  significant,  both  in  time  and  cost. 

ICAPS  water  masses  are  representative  of  very  large  areas 
of  deep  water.  For  application  to  coastal  hydrography  it  is 
better  to  create  a  file  of  the  shelf  area  exclusively.  A 
ten  year  period  may  be  representative,  but  a  longer  period, 
including  as  many  observations  as  possible  should  be  better. 
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The  use  of  a  single  salinity  value  derived  from  a  large 
number  of  observations  is  the  easiest  method  to  be  applied 
in  the  computation  of  depth  corrections.  One  further  step 
is  to  make  a  surface  salinity  observation  and  use  that  value 
for  the  computation  of  the  first  layer  correction.  The  im¬ 
provement  of  this  method  over  the  method  currently  in  use, 
from  a  practical  point  of  view,  is  remarkable. 

This  conclusion  is  applicable  to  the  entire  West  Coast 
of  the  United  States,  excluding  the  Strait  of  Juan  De  Fuca 
and  Puget  Sound  where  the  results  were  not  favorable,  and 
to  San  Francisco  Bay  area  where  there  were  insufficient  data. 

In  summary,  the  applicability  of  the  use  of  historical 
data  in  the  computation  of  depth  corrections  is  as  follows: 

1 .  Northern  Area 

A  single  value  of  salinity  of  31  °/o0  is  applicable 
throughout  the  entire  water  column  except  for  the  surface, 
where  the  salinity  must  be  measured.  In  these  situations, 
the  more  complex  method  4,  wherein  the  historical  salinity 
profile  derived  from  the  ICAPS'  data  file  is  adjusted  to  the 
surface  value,  may  also  be  used. 

2 .  Southern  Area 

The  use  of  a  single  salinity  value  of  33  %0  ,  appli¬ 
cable  throughout  the  entire  water  column,  is  sufficiently 
accurate  for  this  region.  Alternatively,  the  historical 
salinity  curve  derived  from  ICAPS  can  also  be  used  with  no 
adjustment  for  the  surface  layer. 


\ 
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3 .  Strait  of  Juan  De  Fuca  and  Puget  Sound 

The  use  of  historical  salinity  data  derived  from 
of  the  methods  above  is  not  recommended  in  this  area. 
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APPENDIX  A 


WILSON'S  EQUATION 

Wilson  first  presented  his  empirical  formula  for  the 
computation  of  the  speed  of  sound  in  sea-water  in  the  June 
1960  issue  of  the  Journal  of  the  Acoustical  Society  of  .Amer¬ 
ica  [Wilson,  1960a].  In  the  October  issue  of  the  same  year, 
a  slightly  modified  version  of  this  formula  was  introduced, 
with  an  extended  range  of  applicability.  This  version  is 
the  formula  commonly  used  and  accepted  as  best  representing 
the  sound  speed  dependence  on  pressure,  temperature  and 
salinity.  It  fits  data  within  the  ranges  of  -4°C  to  30°C 

-»  .  9 

in  temperature,  1  kg/cm"  to  1 000  kg/cm“  in  pressure  and  0  %o 
to  37  °/0o  in  salinity,  with  a  standard  deviation  of  about 
0.30  m/s  [Wilson,  1960b]4. 

A  value  of  1449.14  is  given  for  the  sound  speed  at  T =  0°C. 

2 

P=  0.0  kg/cm  and  S =  35  %o  •  Corrections  are  added  to  this 
reference  speed  in  order  to  compensate  for  the  variations  in 
pressure,  temperature  and  salinity,  either  alone  or  in  a 
combined  effect.  The  revised  formula  is  as  follows: 


An  analysis  of  this  formula  for  values  of  pressure, 
temperature  and  salinity  outside  the  range  for  which  the 
equation  was  derived  showed  that  the  error  can  go  up  to 
2.13  m/s  for  extreme  conditions  of  temperature  and  salinity 
expected  to  occur  in  the  sea  [Wilson,  1962]. 
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V  =  1449.14  +  VT  +  Vp  ♦  Vs  +  VSTP 


where  = 


V 


P 


V 


S 


V 


STP 


4.5721 T  -  4.4532  x  10~2T2  -  2.6045  x  10'4  T3 
+  7.9851  x  10*6  T4 

1.60272  x  10'1  P  +  1.0268  x  10'5  P2 
+  3.5216  x  10‘9  P3  -  3.3603  x  10‘12  P4 
1. 39799  (S  -  35)  +  1.69202  x  lO'^S-SS)2 
(S  -  35)  (-1.1244  x  10'2  T  +  7.7711  x  10'7  T2 

+  7.7016  x  10‘5  P  -  1.2943  x  10*7  P2 

+  3.1580  x  10~8  PT  +  1.5790  x  10'9  PT2) 

+  P(-1.8607  x  10‘4  T  +  7.4812  x  10'6  T2 
+  4.5283  x  10'8  T3)  +  P2(-2.5294  x  10'7  T 

+  1.8565  x  10'9  T2)  +  P3(-1.9646  x  10'10T) 


The  units  of  temperature,  pressure,  salinity  and  sound 
speed  are,  respectively,  degrees  Celsius,  kilograms  per  square 
centimeter,  parts  per  thousand  and  meters  per  second.  While 
temperature  and  salinity  are  usually  directly  measured, 
pressure  is  computed  from  depth.  The  dependence  of  pressure 
on  depth  is  not  linear  and  several  expressions  may  be  applied 
for  this  computation.  An  iterative  process  is  used  by  NOS 
to  determine  pressure  [Wallace,  1971].  This  method  was 
extracted  from  the  User's  Guide  to  NODC's  Data  Services 
[NOAA,  Department  of  Commerce,  1974],  where  it  is  described 
as  follows: 
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a  =  -  9.5445865  x  10~2  +  8.14876577  x  lO^S 
o 

-  4.8249614  x  10'4  S2  +  6.7678614  x  10‘6  S3 
go  =  0.980616  -  2.5928  x  1O'3(cos20)  +  6.9  x  1O'6(cos20)2 
P(surface)  =  10.1525  decibars 

gi  =  8o  +1-101x  10’7Cdi  +  di-i) 

Pi  -  (1  +  10'3  aT)/R 

R  =  1  -  [4.886  x  10‘6P/(1  +  1.85  x  10~5P)] 

+  P[-2. 2072  x  10'7  +  5.675  x  10'8T  -  6.63  x  10'10T2 
+  4.0  x  10'12T3  +  o  (1. 725  x  10*8  -  3.28xl0~10T 
+  4.0  x  10'12T2)  +  °02('4-5  x  10’ 11  +  1 0  ‘ 12  T3  ] 

+  P2[-6.68  x  10' 14  -  1.24064  x  10~12I  +  2.14  x  10'14T2 
+  oQ(-4.248  x  10'13  +  1.206  x  10_14T  -  2 . 0  x  10'16T2) 

+  oq2  (1.8  x  10' 15  -  6.0  x  10' 17 T)  ]  +  P*(1.5  x  10'17T) 


1. 


2. 


A  first  approximation  gives  PV  and  p^: 

pi  =  +  Pi-igi(di  '  di-i) 

with  PV  computed  this  way,  p/  is  calculated. 

A  second  approximation  using  pV  gives  PV  and  p^ 

pi  ■  pi-i  *  *  pi’ >*itdi  -  di-i> 


with  PV 

l 


p^  is  calculated. 
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5.  Finally,  with  P ^ 

Pi  =  Pi-1  +  1/2  tpi-l  +  pi}gi(di  "  di-l} 
then,  P  *  0.10197  ?i 

In  the  above  expressions,  the  subscript  indicates  the  depth 
that  is  being  computed,  and  the  number  of  primes  is  inversely 
proportional  to  the  accuracy  of  the  approximation.  The  mean¬ 
ing  of  the  different  parameters  is  the  following: 
d  =  depth  in  meters 
T  =  temperature  in  degrees  Celsius 
S  =  salinity  in  parts  per  thousand 
p  =  density 

P  =  pressure  in  kilograms  per  square  centimeter 
0  =  latitude  in  degrees. 
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Figure  B-2.  Location  of  stations  for  the  northern  area 
and  Strait  of  Juan  De  Fuca/Puget  Sound,  in  spring. 
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APPENDIX  C 
COMPUTER  PROGRAM 

PROGRAM  SVC 

PROGRAM  SVC  USES  HISTORICAL^  OERIVED  SALINITY  IN  THE 
COMPUTATION  OF  DEPTH  CORRECTIONS  DUE  TO  VARIATIONS 
IN  SOJND  SPEED.  IT  USES  FOUR  DIFFERENT  METHODS 
FOR  DERIVING  THE  HISTORICAL  VALUES.  THE  DIFFERENT 
DEPTH  CORRECTIONS  ARE  COMPARED  TO  THE  CORRECTION 
08TAINE0  USING  THE  OBSERVED  SALINITY. 


INTEGER  ON,DD,YY»TNQ,TNSO 
REAL  IDL.Ll.Lf ID.LC 

DIMENSION  AVE(43,8 ),AVE1<43  I , COR < 30  )  ,CORl (  30)  ,C0R2(30) 

1  ,C0R3( 30) ,C0R4( 30) , DC0R1 l 3 0 ) , DC0R2 ( 30 ) , DC0R3(30 ) , 

2  DC0R4I  30) ,  DC0RR1  30)  ,  OCR  IK  2000)  ,0CR22(  2000  »  , 

3  DCR3 3 (2000), DC R44< 2 000 ) , ID ( 59 ) , I D1 1  43 ) f L ( 59) ,L I ( 43 ) 

4  , PCI  30) , PC  1(30) , PCI L (2000)  , PC2 ( 30 > , PC22 ( 2000 ) , 

5  PC3(30)»PC33(2000)»PC4(30) ,PC44<2000) ,S(3 0), 

6  SH(22, 8) f SHI l (30) »SH22(3Q) , SH3 { 30  ) , SH4( 30 ) , SHH( 22 ) , 

7  SHHH (2  2) ,T(30),V<30)  ,Vi(30) ,V2(30) ,V3<30) ,V4(30) , 

8  H  22) till  30) 

100  FORMAT  (F7.I.8F7.3) 

101  FORMAT  (312,1 5,AI, I  6, A  1 , IX, A4, I  3  ) 


T  (//LIOX.'X  X'/1UX,' 

1  *X  X'/llOX.^X  X'  ) 

215  FORMAT  ( • 1 •  ,45X, 2 7 (  • +  • ) /46X , • 

1  STATISTICS 

2  46X, 27 ( ' +*  )//  I 

216  FORMAT  ( /  /49X,' PC  1 '  »4X » • PC2 • , 

217  FORMAT  (/•  TOTAL  NUMBER  OF  OS 

1  ,'  NUMBER  OF  OBSERVATIONS  OFF :  •  ,  1 IX , 41 7 / 

2  '  PERCENTAGE:' ,34X,4F7. 2//'  NUMBER  OF  STATIONS: 

3  120/'  NUMBER  OF  STATIONS  OFF :•  ,2  IX, 41 7/ 

4  '  PERCENTAGE:' ,34X,4F7. 2/) 

218  FORMAT  (/'  LARGEST  SINGLE  VALUE  OBSERVED : '  , 1 5 X, 4F 7. 3  / 

1  '  LARGEST  MEAN  VALUE  OBSERVED: ',  17X, 4F7. 3/ 

2  '  LARGEST  STANDARD  DEVIATION  OBS ERVED: ' . 9 X ,4F7 . 3 ) 

219  FORMAT  (//•  MEAN  OF  THE  MEANS  : '  , 2 7X ,4F7. 3/ 1  STANDARD' 

1  ,'  DEVIATION  OF  THE  MEANS: ' t L3X,4F7.4) 


' P  C4  •  » 

.III/'  TOTAL' 
417/ 
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nnoon  on 


FAC (VII =( VI-VCJ/VC 

PCTG(M,N)  =FLOAT(M) /FLOAT (N) *100. 

VC  IS  THE  S0UN0  SPEED  ASSUMED  IN  THE  FATHOMETER. 
VC= 1463 . 04 


I.  COMPUTE  HISTORICAL  SALINITY  AVERAGES 


C 

c 


10 


11 

12 
13 


N=2  2 

DO  10  1=1  ,N 

READ  (5,1001  U  l),(Smi,Ki  ,K=1,8> 

CONTINUE 
WRITE  (6,200) 

WRITE  (6,201)  (Z< I)  ,(SH(I,K» ,K=1,8)  ,1  =  1, N) 

J=0 

M=2*N- 1 
DO  13  K= 1 , 8 
DO  11  1  =  1,  N 

SHH  (  I )  =  SH(  I  ,  K  ) 

CONTINUE 

CALL  MEANH  J ,N , M , Z , SHH , AVE 1 ,  IDl.Ll) 

DO  12  1=1, M 

AVE ( I ,K)  =  AVE1(  I) 

CONTINUE 

J=l 

CONTINUE 
WRITE  (6,202) 

WRITE  (6,203)  ( ID  1 ( I ) , ( AVE( I , K ) , K= 1 , 8 ) , I = 1 , M ) 


0N=0 

TN0=0 

NOF  F1  =  0 

NOF  F2=0 

N0FF3=0 

NOF  F4=0 

DCOR1L=0. 

DC0R2L=0. 

DC0R3L=0. 

DCOR4L=0. 

PC1L=0. 

PC2  L=0. 

PC3L=0. 

PC4L=0. 

dcriml=o. 

QCR2ML=0. 

DCR3ML=0. 

DCR4ML=0. 

PC1ML=0. 

PC2ML=0. 

PC3ML=0. 

PC4ML=0. 

OCR 1SL=0. 

DCR2SL=0. 

0CR3SL=0. 

DCR4SL=0. 

PC  1 SL=0 . 

PC2SL=0. 

PC3SL=0. 

P C4SL=0. 
NSOFF 1=0 
NSOFF2=0 
NSQFF3=0 
NS0FF4=0 
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2.  REAO  DATA 


C 


20 


21 

22 

23 


24 

25 

26 


27 

28 


READ  (5, 101J 
IF  (DD.E3.0) 


WRITE  (6,  204) 
WRITE  (6,205) 
NOFF11=0 
NOF  F22=0 
N0FF33=0 
N0FF44=0 
DCR1TL-0 . 

DCR  2TL-0. 

DCR  3TL=0 . 
DCR4TL=0. 

PC  1 TL  =  0. 

PC2TL=0. 

PC3TL=0. 

PC4TL=0. 

READ  (5,102)  , 

IF  (KONT.EQ.O) 

DO  22  1=2,10 

IF  (Zl(I).NE.O 
N = I  —  1 


YY , MM , DD, LAT »H1,L0NG*H2»NS»K0NT 
GO  TO  80 
NS,  MM,DD ,  YY 
LAT, Hi, LONG, H2 


(Zl( I  )  ,1  =  1,10) 

0  TO  23 


)  GO  TO  21 


GO  TO  26 

IF  (I.EQ.10)  N=I 
CONTINUE 
GO  TO  26 

READ  (5,102  )  (Zl(  I  ), 1  =  11,20) 
DO  25  1=11,20 

IF  (Zl(I).NE.O.) 

N=  I  - 1 
GO  TO  26 
IF  (I.EQ.20)  N= I 
CONTINUE 

READ  (5,102)  (T( I) , 

REAO  (5,102 )  (S(  I ) , 


GO  TO  24 


1  =  1  ,  N) 
1=1, N) 


ON=QN+ 1 

IF  (LAT. GT. 4000)  GO  TO  27 
K=1 

IF  (MM.GE.4)  K=2 
IF  (MM.GE.7)  <=3 
IF  (MM.GE.IO)  K  =  4 
GO  TO  28 
CONTINUE 
K=5 

IF  (MM.GE.4)  K  =  6 
IF  (MM.GE.7)  K  =  7 
IF  (MM.GE.IO)  K  =  8 
CONTINUE 


3.  MAKE  COMPUTATIONS  WITH  TRUE  VALUES 


CALL  LAYER  (N,Z1,L,  ID) 

CALL  CORR  EC ( L  AT , VC , N, 71 ,T , S , L , V, COR ) 


4.  MAKE  COMPUTATIONS  WITH  200,  1000,3000  M  AVERAGES 


SH1=AVE(  19, K) 

IF  (Z1(N) .GT.200.)  SHl=AVEl 33 ,K ) 
IF  (ZUNI  .GT.  1000.  )  SH1  =  AVE(43,K) 
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00  40  1=1 • N 
SH1KI  J=SHl 
CONTINUE 

CALL  CQRRECILAT, VC , N, Z t ,T , SHI  1 , L , V 1 , COR  1 ) 

CALL  OCRR(N»COR1»COR,IO»OCOR1 ,PC1) 

00  41  1=1, N 

AP1  =  AB  S ( PC  1( II ) 

IF  (AP 1.GT.PCITU  PC  1TL= AP 1 
IF  (AP1.GT..25)  NOFF 11=N0F F 1 1  + 1 
CONTINUE 

IF  (PC1TL  .GT.PC1U  PC  IL=PC  ITL 
CALL  STAT  ( N, PC l, PC 1M, PC1S) 

IF  (ABS(PCIM) .GT.PC1ML)  PC1ML=A8S ( PC IM ) 

IF  (ABS(PCIS) .GT.PCtSL )  PC1SL=ABS ( PC1S) 

PC  1 1 (ON )  =  PC  IM 
PRC11=PCTG(N0FF11 ,N  I 
IF  (NOFFll.GT.O)  NSOFF 1=NS0FF 1+ 1 

5.  MAKE  COMPUTATIONS  WITH  50,  100,  150,  200,  300,  5 00, 
750,  1000,  1500,  3000  M  AVERAGES 


SH2=AVE(9,K) 

IF  IZKN)  .GT.50.)  SH2=A VE ( 1 3,K) 

IF  (ZKN)  .GT.  100.  I  SH2  =  AVE(  17, K) 

IF  (ZKN)  .GT.  150.  )  SH2  =  AVE(  19  ,K  ) 

IF  (Z1(N> .GT. 200. )  SH2=AVE( 23 ,K ) 

IF  (ZKN)  .GT.300.  )  SH2=AVE(  27  ,K  ) 

IF  (Zl(N) .GT. 500.)  SH2=AVE( 32 ,K) 

IF  (ZKN)  .GT.750.  )  SH2=AVE(  33  ,K ) 

IF  (ZKN)  .GT. 1000.  )  SH2=AVE  (  3  7  ,K  ) 

IF  (ZKN)  .GT.  1500.  )  SH2=AVE  (  43  ,K  ) 

DO  50  1=1, N 
SH22 ( I ) =SH2 
CONTINUE 

CALL  CORREC(LAT»VC,N,Zi,T»SH22»L,V2»COR2) 

WRITE  (6,206)  SH1,SH2 

CALL  DCRR (N , COR2 , COR, I  0 ,DCOR2 , PC  2 ) 

DO  51  I  =1  »N 

AP2=ABS(PC2( I) ) 

IF  (AP2.GT.PC2TL)  PC2TL=AP2 
IF  (AP2.GT..25)  NQFF 22=N0F F2 2+ l 
CONTINUE 

IF  (PC2TL .GT.PC2L )  PC2L=PC2TL 

CALL  STAT  (N, PC2,PC2M, PC2S) 

IF  (ABS(PC2M) .GT.PC2ML)  PC2ML= ABS ( PC2M ) 

IF  (ABS(PC2S) .GT.PC2SL)  PC2 SL=A BS ( PC2S ) 

PC22(ON)=PC2M 

PRC  22=PCTG (NOFF2  2  ,  N ) 

IF  ( NOFF22.GT .0)  NSOFF2=NSQFF2+l 


MAKE  COMPUTATIONS 
SALINITY  VALUE 


EACH  CORRESPONDING  HISTORICAL 


DO  60  1=1,22 

SHHH(  I )  =  SH  ( I ,  K  ) 

CONTINUE 

CALL  INTRPL  ( 22, Z , SHHH,N,Z1 ,SH3 ) 

CALL  C0RREC(LAT,VC,N,Z1 ,T , SH3 , L , V3 , C0R3 ) 
CALL  DC RR(N,C 0R3, COR, ID, DC0R3»PC3 ) 

00  61  1=1 ,N 

AP3= AB  $ ( PC  3( I)  ) 

IF  ( AP3.GT.PC3TL)  PC3TL=AP3 
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IF  (AP3.GT..25) 
CONTINUE 

IF  (PC3TL.GT.PC3L) 


N0FF33=N0FF33+1 


IF  (PC3TL.GT.PC3L)  PC3L=PC3TL 

CALC  STAT  (N»PC3»PC3M,PC3SI 

IF  { ABS  ( PC3M) .GT . PC3ML )  PC3ML = ABS ( PC3M ) 

IF  ( ABS  (  PC3  S) .GT.PC3SL)  PC3 SL=ABS ( PC3S ) 

PC33 (QN)=PC3M 

PRC33=PCTG ( NOFF33  ,  N ) 

IF  (N0FF33.GT .0)  NS OFF 3=NS0FF  3+1 


7.  MAKE  COMPUTATIONS  WITH  ADJUSTED  HISTORICAL 
SALINITY  CURVE 

ADJ  =  S(  1J-SH3I  II 
DF=30.0-Zll  U 
DO  71  1*1 »N 

IF  (ZKII.GE.30.01  GO  TO  70 
SH4(I)=SH3(I)+ADJ*(30.0-Z1(I I  I /DF 
GO  TO  71 

I  SH4( I I =SH3 ( I ) 

.  CONTINUE 

CALL  C0RREC(LAT,VC,N,Z1 , T ,SH4 ,L , V4 ,  CQR4 ) 

CALL  DCRR (N,C0R4, COR, I D,DC0R4,PC4) 

00  72  1=1 »N 

AP4=ABS(PC4( I) I 
IF  (  AP4.GT.PC4TL)  PC4TL=AP4 
IF  (AP4.GT..25)  N0FF44=N0FF44+ l 
>  CONTINUE 

IF  (PC4TL.GT.PC4L I  PC4L=PC4TL 
CALL  STAT  ( N  ,  PC4,  PC4M,  PC4S) 

IF  (ABSIPC4M) .GT.PC4MLJ  PC4ML=ABS( PC4MI 
IF  ( ABS ( PC4  S )  .GT.PC4SL1  PC4SL  =  ABS ( PC4SI 
PC 44 (ON) =  PC4M 
PRC44=PCTG(N0FF44,N) 

IF  (N0FF44.GT .0 )  NS0FF4*NS0FF4+L 
WRITE  (6,207) 

WRITE  (6.208)  ( Z 1 (  I  ) , I D ( I ) , T(  I  )  , S  ( I) ,  SH3(  I ) , SH4( I ) , 

1  V< I ) ,Vl (I ) ,V2( I) *  V3 1  I ) , V4( I  ) ,  I  =  l,N) 

WRITE  (6,209) 

WRITE  (6,210)  ( Z1 ( I ), I D ( I),CORI I ) ,C0R1(I) ,C0R2( I) , 

1  C0R3  ( I  I  ,CQR4(  I  )  ,  DCORK  I)  ,DC0R2(  I  ),DC0R3(  I )  ,DC0R4(I  ) 

2  ,PC1(I ),PC2( I),PC3( I),PC4( I),I=1,N) 

WRITE  (6,211)  N 

WRITE  (6,212)  NOFF 11, N0FF22  »NQFF33  »  N0FF44 , PRC  1 1 , 

1  PRC22,  PRC33,=>RC44 

WRITE  (6,213)  PC1TL,PC2TL,PC3TL,PC4TL,PC1M,PC2M, 

1  PC3M.PC4M, PC1S ,PC2S,PC3S,PC4S 
NOFF=NOFF  1H-NOFF22  +  NOFF33+NOFF44 
IF  (NOFF.NE.O)  WRITE  (6,214) 

TNO=TNO+N 

NQFFI=N0FF1+N0FFI 1 
NOFF2=NOFF2+NOFF22 
NOFF3=NOFF3+NOFF33 
NQFF4=NQFF4+N0FF44 
GO  TO  20 


COMPUTE  STATISTICS 


80  CONTINUE 

PRC l=PCTG(NQFFl , TNO I 
PRC2=PCTG (N0FF2, TNG  ) 
PRC3*PCTG (N0FF3 ,TNO ) 
P  RC4*PCTG (N0FF4, TNO  ) 
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PRCT1=PCTG(NS0FF1,0N) 

PRC T2=f>CTGiNSOFF2f  ON) 

PRC  T 3= PCI G(NSQFF3»ON) 

PRCT4=PCTG ( NS0FF4 ,  ON  I 

CALL  STAT  < ON  ,PC 1 1 , PC  1 1 M, PC H $ > 

CALL  STAT  ( ON  ,PC22 , PC22M, PC  22 S J 
CALL  STAT  l ON  ,PC33 , PC33 M, PC 33 S ) 

CALL  STAT  ( ON  ,PC44, PC44M, PC 44 S) 

WRITE  (6,215) 

WRITE  (6,216) 

WRITE  (6,217)  TNO , NOFF  1  ,NOF F2 , NOF F3,N0FF4,PRC1,PRC2, 
1  PRC3,PRC4,QN,NS0FF1,NS0FF2 ,NS0FF3,NS0FF4, PRCT1, 

0  DP  r  T  ?  .  DP  C  T  3  .  P  B  f  Ti 

WRITE  (6,218)  ^C 1L , PC2L  »PC3L,PC4L,PC 1ML , PC2ML  , 

1  PC3ML  »  PC4ML, PC  1 SL , P  C  2SL , PC3SL , P  C4SL 

WRITE  (6,219)  PC11M,PC22M,PC33M,PC44M,PC11S,PC22S, 

1  PC33S, PC44S 
STOP 
ENO 
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SUBROUTINE  MEAN! 


SUBROUTINE  MEANI  (  J  ,N, K  , Z  ,  P  ,  A VE  ,  I D  ,L  ) 

SUBROUTINE  MEAN!  MAKES  A  WEIGHTED  AVERAGE  OF  PARAMETER 
Pt  USING  SUBROUTINE  LAYER!  FOR  COMPUTATION  OF 
THE  WEIGHTS. 


MEANING  OF  THE  PARAMETERS: 

INPUT: 

"J"  -  FLAG;  AVOIO  RECQMPJTAT I  ON  OF  THE  WEIGHTS  IF 

THEY  HAVE  BEEN  COMPUTED  BEFORE. 

=0  IF  SUBROUTINE  HAS  NOT  YET  BEEN  USED 
WITH  THE  SAME  WEIGHTS. 

=1  (OR  DIFFERENT  THAN  ZERO)  IF  SAME 
WEIGHTS  ARE  TO  APPLY. 

"N"  -  NUMBER  OF  OBSERVATIONS. 

"K"  -  =2N-L  -  DIMENSION  OF  ARRAYS  AVE,  ID,  L 

"Z"  -  ARRAY  FOR  THE  DEPTH.  DIMENSION  N. 

Z(l)  MUST  BE  ZERO. 

»P"  -  ARRAY  FOR  THE  QUANTITY  TO  BE  AVERAGED. 

DIMENSION  N. 

OUTPUT: 

"AVE”-  AVERAGE  ARRAY.  WILL  CONTAIN  THE  WEIGHTED 
AVERAGE  FOR  EACH  DEPTH  GIVEN  AND  ALSO 
FOR  EACH  INTERMEDIATE  DEPTH.  DIMENSION  K. 
"ID"  -  INTERMEDIATE  DEPTH.  DIMENSION  K  (SEE 
SUBROUTINE  LAYER!). 

"L"  -  LAYERS.  DIMENSION  K  (SEE  SUBROUTINE 

LAYER  1 ) . 

REAL  IO.L 

DIMENSION  Z(N),P(K) , AVE ( K > , ID ( K ) , L( K ) , R( 43 ) 

IF  (J.NE.O)  GO  TO  300 

CALL  LAYER1(N,K,Z,L,I0) 

300  CONTINUE 
SUM=O.DO 
AVE ( 1 ) =P ( 1) 

DO  301  1=2, K 

J=IFIX(FL0AT(I  +  1  )/2.  ) 

R ( I ) =P  ( J  ) 

SUM=S J M+L (  1-1 1 *R(  II 
AVE( I ) =SUM/ID( I ) 

301  CONTINUE 
RETURN 
END 
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SUBROUTINE  LAYERL 


SUBROUTINE  LAYER1  (N»K,Z.L»ID) 

SUBROUTINE  LAYERL  COMPUTES  THE  INTERMEDIATE  DEPTHS 
OF  A  GIVEN  WATER  COLUMN  AND  THE  LAYERS  BETWEEN 
THEM.  IT  OIFFERS  FROM  SUBROUTINE  LAYER  BY  ,tc 
CONSIDERING  EACH  GIVEN  DEPTH  AS  AN  INTERMEDIATE 
OEPTH. 

MEANING  OF  THE  PARAMETERS: 
f  Mp UT  « 

*  i*N»  -  NUMBER  OF  DEPTHS  GIVEN.  , 

11*11  -  =2N-1  -  NUMBER  OF  INTERMEDIATE  DEPTHS. 

DIMENSION  OF  ARRAYS  L»  ID. 

'•Z"  -  ARRAY  FOR  THE  DEPTH.  DIMENSION  N. 

Z(L»  MUST  BE  ZERO. 


OUTPUT: 

t«  - 

tun**  - 


REAL  I D» L 
DIMENSION 
M=N-1 
DO  401  1= 
J=2*I 
I D( J  I  * 
IF  (ID 
CONTINUE 
DO  402  I* 
J=2*I- 
10( JJ= 
CONTINUE 
LAYERS 
L(1)=ID(2 
J  =  J-l 
DO  403  1= 
L  (  I » =  I 
CONTINUE 
RETURN 
END 


LAYERS  IN  BETWEEN  EACH  INTERMEDIATE  DEPTH, 
K-l  LAYERS  ARE  CONSIDERED. 

INTERMEDIATE  DEPTHS.  BOTH  EACH  GIVEN 
DEPTH  AND  HALF-WAY  IN  BETWEEN  TWO 
SUCCESSIVE  GIVEN  DEPTHS  ARE  CONSIDERED. 
ID ( L )  IS  THE  SURFACE  AND  ID( K ) THE  LAiT 
GIVEN  DEPTH.  DIMENSION  K. 

ZCN) ,ID(KI tL(K» 


<z(  n+z(  i+i)  )/2. 

( J ). EQ. 700 . )  I D( Jl *7  50 . 


2*  J 

D(  I  +  U-IDI  I) 
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SUBROUTINE  LAYER 


SUBROUTINE  LAYER  (N,Z,L,IDi 

SUBROUTINE  LAYER  COMPUTES  THE  INTERMEDIATE  DEPTHS 
OF  A  GIVEN  WATER  COLUMN  AND  THE  LAYERS  BETWEEN 
THEM.  IT  DIFFERS  FROM  LAYER1  BY  NOT  CONSIDERING 
THE  GIVEN  DEPTHS  AS  INTERMEDIATE  DEPTHS. 

MEANING  OF  THE  PARAMETERS: 

INPUT: 

"N"  -  NUMBER  OF  DEPTHS  GIVEN. 

«Z"  -  ARRAY  FOR  THE  DEPTH.  DIMENSION  N. 

Z(l)  MUST  BE  ZERO. 


OUTPUT: 

"L"  -  ARRAY  FOR  THE  LAYERS.  DIMENSION  N. 
"ID"  -  ARRAY  FOR  THE  INTERMEDIATE  OEPTHS. 
DIMENSION  N. 

REAL  ID,L 

DIMENSION  ZINI,  L ( N  )  »  ID(N) 

FIND  INTERMEDIATE  DEPTHS  ID 
M=N— 1 

00  501  I  =  L  ,  M 

I  D(  I  )  =  I Z  (  IHI+ZI  1)1/2. 

501  CONTINUE 
ID  (  N I  =Z  I  N  ) 

FIND  LAYERS 


L(1»*I0<1 I 
DO  502  1=2, N 

L(II  =  ID( II-IOI  I- II 
502  CONTINUE 
RETURN 
ENO 
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SUBROUTINE  CORREC 


SUBROUTINE  CORREC  (  LAT  ,  VC  ,N  , Z  , T ,  S » L ,  V, C OR) 

SUBROUTINE  CORREC  COMPUTES  THE  DEPTH  CORRECTIONS 
DUE  TO  VARIATIONS  IN  SOUND  SPEED. IT  USES 
SUBROUTINE  WILSON  TO  COMPUTE  SOUND  VELOCITY 
AT  DEPTH  Z,  ANO  APPLIES  THE  CORRECTIONS  TO  THE 
ASSUMED  SOUND  SPEED  VC.  CORRECTIONS  ARE  DETERMINED 
BY  LAYERS  L  AND  SUMMED  ALGEBRAICALLY  TO  GIVE  THE 
CORRECTION  AT  THE  BOTTOM  OF  EACH  LAYER. 

MEANING  OF  THE  PARAMETERS: 

INPUT: 

"LAT"-  LATITJOE  IN  DEGREES. 

"VC"  -  SOUNO  SPEED  ASSUMED  IN  THE  FATHOMETER. 

"N"  -  NUMBER  OF  DEPTHS  GIVEN. 

"Z"  -  ARRAY  FOR  THE  DEPTH.  DIMENSION  N. 

Zlll  MUST  BE  ZERO. 

"T"  -  ARRAY  FOR  THE  TEMPERATURE.  DIMENSION  N. 

"S"  -  ARRAY  FOR  THE  SALINITY.  DIMENSION  N. 

OUTPUT: 

"V"  -  ARRAY  FOR  THE  SOUND  SPEED.  DIMENSION  N. 

"COR"-  ARRAY  FOR  THE  CORRECTIONS.  DIMENSION  N. 

REAL  L.LC 

DIMENSION  Z(N)«T(N),S(N),L(N),  COR ( N  )  »V( N ) 

FAC (VI ) =( VI-VCJ/VC 

CALL  WILSON  l  LAT .  Z  (  1 ) ,  T  ( l )  ,  SI  l )  ,  V  (  l ) ) 

LC  =  L( l ) *F  AC ( V  (1) ) 

CQR(U=LC 
DO  600  1=2, N 

CALL  WILSON  (LAT.Zm,T(  l)  ,S(I),V(I)  ) 

LC=L ( I ) *FAC( V ( I )  I 
COR  < II =COR  <1  —  11 ♦LC 
600  CONTINUE 
RETURN 
END 
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SUBROUTINE  OCRR 


SUBROUTINE  OCRR  l N , COR  1  ,C0R 1 1 D, OCOR , PC ) 

SUBROUTINE  DCRR  GIVES  THE  DIFFERENCES  BETWEEN 

THE  DEPTH  CORRECTIONS  OBTAINED  BY  A  PARTICULAR 
METHOO  AND  A  REFERENCE  METHOD. 

MEANING  OF  THE  PARAMETERS: 

INPUT: 

"N"  -  NUMBER  OF  OEPTHS  GIVEN. 

"COR1"-  ARRAY  FOR  THE  CORRECTIONS  FROM  THE 
METHOD  USED.  OIMENSION  N. 

"COR"  -  ARRAY  FOR  THE  CORRECTIONS  FOR  THE 
REFERENCE  METHOD.  DIMENSION  N. 

"ID"  -  INTERMEDIATE  DETHS  TO  WHICH  THE  CORRECTIONS 
APPLY.  DIMENSION  N. 

OUTPUT: 

"OCOR"-  ARRAY  FOR  THE  DIFFERENCE  IN  THE 
CORRECTIONS.  DIMENSION  N. 

"PC"  -  PERCENTEGE  OF  THE  DIFFERENCE  TO 
THE  DEPTH. 

"PC"  -  ARRAY  FOR  THE  PERCENTAGE  OF  THE 

DIFFERENCE  TO  THE  QEPTH.DI MENSION  N. 

REAL  ID 

DIMENSION  CORUN)  ,  COR(  N  )  t  ID  { N  ) ,  DCOR  {  N) ,  PC  (  N  ) 

DO  700  I*1,N 

OCORU  )*C0R1(  I)-COR(  I) 

pcu)=dcor  m/io<  n*ioo. 

700  CONTINUE 
RETURN 
END 
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SUBROUT I  ME  STAT 


SUBROUTINE  STAT  (  M , A, ME  AN, STQEVJ 

SUBROUTINE  STAT  COMPUTES  THE  MEAN  AND  STANDARD 
DEVIATION  OF  THE  ARRAY  A,  WHICH  HAS  M  ELEMENTS 

DIMENSION  MUST  BE  ASSIGNED  IN  MAIN  PROGRAM  FOR  ARRAY  A 

DIMENSION  A (M I 

REAL  MEAN 

COMPUTE  MEAN 

IF  (M.GT.1J  GO  TO  800 
MEAN=A< i> 

STDEV=0 . 

RETURN 

800  SUM=0. 

□0  801  1=1, M 
SUM=SUM*A(  I) 

801  CONTINUE 

MEAN=SUM/ FLOAT!  M  J 

COMPUTE  STANDARD  DEVIATION 

SUM=0. 

DO  802  J=  1 ,  M 

TERM*!  A( J I -MEAN  1**2 
SUM=SJM+TERM 

802  CONTINUE 

STDEV=SQR  T ( SUM/FLOATI M-l ) ) 

RETURN 

END 
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SUBROUTINE  WILSON 


SUBROUTINE  WILSON  ( LAT  ,  Z , T, S , C ) 

SUBROUTINE  WILSON  COMPUTES  THE  SOUND  VELOCITY  IN  WATER 
USING  WILSON'S  EQUATION  (I960).  THE  EQUATIONS  USED 
IN  THIS  SUBROUTINE  WERE  EXTRACTED  FROM  "USER'S 
GUIDE  TO  NODC'S  DATA  SERVICES  ( NOAA , OE P ARTMENT 
CF  COMMERCE, 1974). 

THE  PRESSURE  USED  IN  THIS  EQUATION  IS  COMPUTED  IN 
AN  ITERATIVE  PROCESS  ALONG  DEPTH.  THE  SUBROUTINE 
SHOULD  BE  CALLEO  SEQUENTIALLY  FOR  THE  DIFFERENT 
DEPTHS,  STARTING  FROM  THE  SURFACE  C  Z  =  0  >  - 

MEANING  OF  THE  PARAMETERS: 

INPUT: 

"LAT"-  LATITUDE  IN  DEGREES. 

**Z"  -  DEPTH. 

"T"  -  TEMPERATURE. 

"S"  -  SALINITY. 


OUTPUT: 

"C" 


-  SOUND  SPEED. 


T2=T**2 

T3=T**3 

T4=T**4 

S2=S**2 

S3=S**3 

IF  (Z.NE.O.)  GO  TO  900 
PI=10. 1325 
P=P1 

0LAT=FLQAT(LAT)/100. 

X=COS( 2. *0LAT*3 .1415927/180. ) 

Zl=0. 

900  CONTINUE 

COMPUTATION  OF  DENSITY 

F l=-( T-3. 98>**2*<T+283.)/t503.57*(T*67.25) J 

F2= I.0843E-6*T3-9. 8 18 5E-5*T2* 4. 7867E-3*T 

F3  =  l.667E-8*T3-8.  I6  4E-7*T2  +  1. 803E-5*T 

F= 6. 76786 I36E-6*S3- 4. 82496l40E-4#S2+8. 14876577E-l*S 

DC  1*3. 89 541 4E -2 

DC2=-0.  22  5845  86 

SGMT «Fl ♦( F+DC1)*( 1 «-F2+F3*( F+DC2 )  ) 

COMPUTATION  OF  PRESSURE 

SGM0=-9.3445863E-2+  8.14876577E-I*S-4.8249614E-4*S2+ 
1  6.76786L4E-6*S3 

SSMQ2=SGM0**2 

G0=0.  980616-2.592  3E-3*X'-6.9E-6*X**2 
GI=G0+i.l01E-7*(Z*Zl) 

I  a0 
TP  =  Pl 

901  TP2-TP**2 
TP3=TP**3 

R*i.-(4.886E-6*TP/(  l.  +  l  .83E-5*TP  I  )+TP*(-2.2072E-7* 
3.673E-8*T-6.53E-10*T2+4.E-12*T3+SGM0*<  1.7  25E-8- 


1 

2 
3 


3.28E-  10*T*4.0E-12*T2)  +SGM02*  (-4.5E-11+  L.  E-l 2*T)  )  + 
TP2*(-6.b8E-14-l  .24064E-12*T+2. 14E-  14* T2  + 
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4  SGM0*( -4.248E-13+l.206E-14*T-2.E-16*T2)+ 

5  SGM02*(1.8E-15-6.0E-17*T) > +  TP3 *1 . 5E- 17*T 
R3= 1 1 - > 1 - E-3*  SGMT ) /R 

IF  (Z.EQ.O.I  GO  TO  903 
1  =  1  +  1 

IF  (I.NE.  L)  GO  TO  902 
TP=P 1+RO 1*GI* ( Z-Z 1 ) 

GO  TO  901 

902  T°=P1+.5*(R01+R3)*GI*(Z-Z1) 

IF  (I .EQ«  2 )  GO  TO  901 

903  P=. 10197*TP 

SOUND  SPEED  EQUATION 


C 

c 


P2=P**2 
P3=P**3 
P4=  P**4 

VP  = 1.60272E-1 *P+i . 0263E-5*P2+3. 5216E-9*P3 
1  -3.3603E-12+P4 

V$=l. 39799* JS-35. ) + 1 . 69202E-3 *( S-3  5.  )**2 

VT  =  4. 5721 *T-4.453 2E-2* T2-2. 6045E-4*T3  +  7. 98  51 E-6*T4 

VSTP=(S-3  5. »*(-i.  1244E- 2*T+ 7.  77 1 1 E- 7*T2  + 7 . 70 1 6E-5*P 

1  -l.2943E-7*P2+3 . 1 53 E-8 *P *T +1 . 57  9E-9*P*T2  ) 

2  +P*(- t.8607E-4*T+7. 4312E-6*T2+4. 5233E-8*T3 ) 

3  +P2*<-2. 5294E-7*T+l .8563E-9*T2 ) +P3*( - 1 . 9646E-1 0* T ) 

C= 1449. 14+VP+VS+VT+VSTP 


R01=R0 

RETURN 

END 


8  0 


ooononnnnonnonnnnnnnonnnnno  <  >oor.o 


SUBROUTINE  INTRPL 


SUBROUTINE  I NTRPL ( L , X  ,Y , N, U, V ) 

SUBROUTINE  INTRPL  IS  PART  OF  THE  IBM  SCIENTIFIC 
SUBROUTINE  PACKAGE  EXISTING  AT  THE  NAVAL 
POSTGRADUATE  SCHOOL  COMPUTER  CENTER. 

INTERPOLATION  OF  A  SINGLE-VALUED  FUNCTION 

THIS  SUBROUTINE  INTERPOLATES,  FROM  VALUES  OF  THE  FUNCTION 
GIVEN  AS  ORDINATES  OF  INPUT  DATA  POINTS  IN  AN  X-Y  PLANE 
AND  FOR  A  GIVEN  SET  OF  X  VALUES  (ABSCISSAS),  THE  VALUES  OF 
A  SINGLE-VALUED  FUNCTION  Y  =  Y( X) . 

THE  INPUT  PARAMETERS  ARE 


L  = 


X  = 


Y  = 
N  = 


U  = 


THE 
POINTS 


X  VALUES 


Y  VALUES 


NUMBER  OF  INPUT  DATA  POINTS 
(MUST  BE  2  OR  GREATER) 

ARRAY  OF  DIMENSION  L  STORING 
(ABSCISSAS)  OF  INPUT 
(IN  ASCENDING  ORDER) 

ARRAY  OF  DIMENSION  L  STORING  THE 
(ORDINATES)  OF  INPUT  DATA  POINTS 
NUMBER  OF  POINTS  AT  WHICH  INTERPOLATION 
Y  VALUE  (ORDINATE)  IS  DESIRED 
(MUST  BE  I  OR  GREATER) 

ARRAY  OF  DIMENSION  N  STORING  THE  X  VALUES 
(ABSCISSAS)  OF  DESIRED  POINTS 
OUTPUT  PARAMETER  IS 

V  =  ARRAY  OF  DIMENSION  N  WHERE  THE  INTERPOLATED 
VALUES  (ORDINATES)  ARE  TO  BE  DISPLAYED 
DECLARATION  STATEMENTS 


OF  THE 


THE 


DIMENSION 
EQUIVALENCE 
REAL 

EQUIVALENCE 


PRELIMINARY 
10  LO=L 

LM1 =L0-1 
LM2=LMl-l 
LPl=LO*l 
NO=N 


X  ( L )  ,  Y  (  L )  ,  U  (  N )  ,  V  ( N  ) 

(P0,X3) ,( Q0»Y3) , (Q1,T3> 

MI, M2, M3, M4,M5 

(UK»DX)»(  IMN,X2»A1»M1) »(IMX,X5,A5»M5)» 
( J  ,  SW , SA ) ,(Y2»W2,W4»Q2) ,(Y5,W3,Q3) 
PROCESSING 


IF ( LM2.LT .0) 

GO 

TO 

90 

IF(NO.LE.O) 

DO  11  1=2, LO 

GO 

TO 

91 

IF(X(  I-l)-X(I)  ) 

11 

,95 

,96 

II 

CONTINUE 

IPV=0 

MAIN 

DO-LOOP 

DO  80  K=  1 ,  NO 

UK=U( K) 

ROUTINE  TO  LOCATE  THE 

DESIRED 

POINT 

20 

I F ( LM2 . EQ . 0 ) 

GO 

TO 

27 

I F( UK . GE. X( LO) ) 

GO 

TO 

26 

IF(UK.LT.X{  D) 

I  MN=2 

I MX=LO 

GO 

TO 

25 

21 

I=( IMN+IMX) /2 
IF(UK.GE.X(  I)  ) 

GO 

TO 

23 

22 

I  MX=  I 

GO  TO  24 

23 

I MN=I +1 

24 

I F ( IMX.GT.IMN) 

I  =1  MX 

GO  TO  30 

GO 

TO 

21 
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I  =1 

GO  TO  30 

I  =L?  1 
GO  TO  30 
1=2 

IF  I  =  IPV 
IF( I.E3 .IPV) 

I  PV=  I 

ROUTINES  TO  PICK  UP  NECESSARY  X  AND  Y  VALUES 
TO  ESTIMATE  THEM  IF  NECESSARY 

J  =  I 


2o 
27 

CHECK 

30 


GO  TO  70 


AND 


AO 


4l 


42 

43 
45 


4o 

47 


48 


IFC  J.EQ.U 
IFU.EQ.LPI  ) 

X  3=X ( J- I ) 

Y  3= Y( J- U 
X4=X( J ) 

Y4=Y( J  ) 

A3=XA-X  3 
M3=(  Y4-Y3 ) / A3 
I F ( L  M  2 . EQ.O) 

I  F( J . EQ  .2  ) 

X2=X( J-2) 

Y2= Y(  J- 21 
A2  =  X3-X  2 
M 2  =  (  Y3-Y2) /A 2 
IF(J.Ea.LO> 

X5  =  X(  J+l) 

Y  5  =  Y  (  J+-  1) 

A4=X5-X4 

M4=(Y5-Y4J/A4 

IFU.E3.2) 

GO  TO  45 
M4=M3+M3-M2 
GO  TO  45 
M2=M3 
M  4=M3 

IF( J.LE.3) 

A l  =  X2-X  C  J  —  3 )  #  , 

Ml=(  Y2-YI  J-3)  )/AI 
GO  TO  47 
M1=M2*M2-M3 
IF( J.GE.LM1 > 

A5=X( J+2J-X5  r 
M5=(  Y(  J  +  2)-Y5)  /  A  5 
GO  TO  50 
M5-M4+M4— M3 


J=  2 
J  =  L  0 


GO 

GO 


TO  43 
TO  41 


GO  TO  42 


M2=  M3  +  M3- M4 


GO  TO  46 


GO  TO  48 


C  NUMERICAL  DIFFERENTIATION 


50 


GO  TO  52 


GO  TO  5 1 


51 

52 


53 


IFU.EQ.LPI  ) 

W2=ABS<  M4-M 3) 

W3=ABS(  M2-M1) 

SW=W2+Q3 
I  F  (  SW.NE.O.  0) 

W  2=0 . 5 
W3=0. 5 
SW=1.0 

T3=(W2*M2+W3*M3 I/SW 

IF( I.E3.1)  GO  TO  54 

W3=A8S(  M5-M4) 

W4=ABS ( M3 -M2 ) 

I FISwtNE.O. 0)  GO  TO  53 

W 3=0.  5 
W  4=0.  5 

T4  = ( W3*M3+W4*M4 ) / SW 

I F ( I .NE.LP1 )  GO  TO  60 

T3=T4 
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S  A=A2-*-A3 

T4*0.  5*( M4+M5-A2*( A2- A3 ) * < M2-M3  )/(SA#SA)  ) 

X3=X4 
Y3  =  Y4 
A3=A2 
M3=M4 
GO  TO  6  0 
54  T  4=13 

SA=A3*A4 

T2  =  0.5*(Mi«-M2-A4*  { A3- A4) * ( M3-M4 ) / l  SA*SA  )  ) 

X3-X3-A4 

Y3=Y3-M2*A4 

A3=A4 

M3=M2 

C  OCT ERMINATIQN  OF  THE  COEFFICIENTS 
60  Q2*12.0*(M3-T3I *M3-T4)/A3 

Q3=(-M3-M3*T3+T4>/( A3*A3) 

C  COMPUTATION  OF  THE  POLYNOMIAL 
70  QX=UK-PO 

80  V(K)=QO+OX*(Q1*QX*(Q2+OX*Q3) ) 

RETURN 
C  ERROR  EXIT 

90  WRITE  (6,2090) 

GO  TO  99 

91  WRITE  (6,2091) 

GO  TO  99 

95  WRITE  (6,2095) 

GO  TO  97 

96  WRITE  (6,2096) 

97  WRITE  (6,  2097)  I  ,  X(  I) 

99  WRITE  (6,2099)  LO,NO 

RETURN 

C  FORMAT  STATEMENTS 

2090  FORMAT l IX /22H  ***  L  «  1  OR  LESS./) 

2091  FORMAT ( IX /22H  ***  N  *  0  OR  LESS./) 

2095  FORMAT ( 1X/27H  ***  IDENTICAL  X  VALUES./) 

2096  FORMAT ( 1X/33H  ***  X  VALUES  OUT  OF  SEQUENCE./) 

2097  FORMAT { 6H  I  *, I  7 , 10X, 6HX( I)  =,E12.3> 

2099  FORMAT ( 64  L  =,I7,L0X,3HN  =,17/ 

L  36H  ERROR  DETECTED  IN  ROUTINE  INTRPL) 

END 
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